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' Abstract: Resummed results are discussed for DIS event shape variable distributions 

' involving the thrust defined using the thrust axis of the current hemisphere, the C para- 

■ meter and the jet mass. The resummation of these variables requires the development of 

CO i some techniques different to those applied before for the resummation of DIS event shapes. ^ ^ 

Including power corrections and matching a comparison with HI data is presented. 



1. Introduction 



^ The study of event shape variable distributions is a rather well developed area within 

', perturbative QCD. It was observed nearly a decade ago (see e.g. ||]) that shape variable 

I I distributions, although infrared and collinear safe observables suffer from the presence of 

logarithms that spoil the convergence of the fixed order expansion. For a typical event 
' shape cross section, which is the integral of the distribution over a limited range of the 

shape variable V, the leading behaviour at n^^ order in perturbation theory is 

ii(y)~a:in2"i + ... (1.1) 

The dots represent less singular terms which, however, are still relevant for a description of 
the data in the kinematic region V <^ 1, the two-jet (1+1 jet ) limit in e~^e~ annihilation 
(DIS). These logarithms can, in many cases, be resummed into an expression that is both 
meaningful (physically) and is practically useful for comparisons with data. The state of 
the art resummation is generally till single logarithmic accuracy, i.e. aims to account for 
all terms more singular than and upto a" In" y . 

Resumming the double logarithms (Eq. ^?^) is straightforward and involves the imple- 
mentation of a probabilistic parton branching evolution pattern on each jet which method 
also treats collinear enhancements which are sources of single logarithms. However, in 
certain cases the resummation of single logarithms can be more involved. Recent examples 
are provided by three-jet observables in e'^e~ annihilation [|| and by a class of observables 
sensitive to radiation in only a portion of phase space |^]. The observables we discuss here 
are of this latter kind. 
'Speaker. 
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2. Definitions 



Next we shall define the various event shapes that have most recently been resummed. 
These are the thrust wrt the thrust axis in the Breit frame, the jet mass and the C 
parameter, defined respectively as 

\Pi.n\ 

T = max ^^'^ ' I ' , (2.1) 

^^3E.,^J^.||n|sin^...^ (2.3) M 

where Pi are final state hadron momenta and the sums count only particles in Hq^ the CD 
current hemisphere. Although the above definitions are in terms of the hadron momenta, | — i 
perturbative calculations are performed by replacing the hadron momenta by parton mo- 
menta. This procedure costs us non-perturbative effects that vary inversely as the hard 
scale Q of the process (power corrections) and which we estimate through renormalon based 
techniques. Another definition of the thrust involves using the photon axis in the Breit 
frame and was resummed earlier Q 

= 2 ^^cl^^-^7l _ ^2.4) 

We neglect Z boson exchange although this could be relevant at higher Q values. The 
resummed results (form factors) will in any case be identical to the photon case which 
means that the dominant logarithmically enhanced terms are the same. The fixed order 
results with which one combines the resummed results could be affected through their non 
logarithmic terms. Unfortunately the commonly used fixed order Monte-Carlo programs 
do not implement Z exchange and as such we are unable to correct for this effect. 



3. DIS event shape kinematics 

In this section we discuss the main kinematic differences between the observables resummed 
in this article and those involving the photon axis in the Breit frame [Q]. 

Consider a hard parton at a small angle 9 to the photon axis {z axis) in the Breit frame. 
This angle is most generally acquired through recoil against emissions entering either hemi- 
sphere (current or remnant). Further consider switching off all secondary radiation in the 
current region. In other words, to start with, just imagine a situation where the emissions 
are confined to the target hemisphere and the only object in the current hemisphere is the 
hard parton. Note that one imposes a minimum energy requirement in the current region. 
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which guarantees that the current region is not empty and the shape variables are infrared 
safe. Then one obtains for the thrust defined with respect to the photon axis 



where kt and E are the parton transverse momentum and energy. On the other hand ah 
the variables defined in the previous section vanish (1 — T = p = C = 0), unless there are 
emissions in Tic- This leads to two effects 

• If one wants only values of 1— near zero, which is the region requiring resummation, 
one notices that the kt generated by space-like emissions in the target region has to 
be small. This is not the case for the other variables here, for example the jet mass 
still vanishes even though one may have a large kt (wrt photon axis) single parton 
in the current hemisphere. Placing a restriction on the space-like evolution means 
that the scale of the structure function is not but rather t-^Q"^ . Hence there are x 
dependent single logarithms that are resummed by DGLAP evolution. For the other 
observables the appropriate scale is still Q^. 

• Even though one directly counts only partons (hadrons) in the current hemisphere in 
the definition of all our observables, the thrust wrt the photon axis and normalised 
to Q/2 is a global observable. It counts soft /collinear secondary gluons in either 
hemisphere on an equal footing because of recoil effects discussed above^. On the 
other hand the current jet-mass, the C parameter and the thrust wrt the actual 
thrust axis are non-global observables sensitive to radiation in only a part of phase 
space, the current hemisphere. This necessitates resumming single logs generated by 
coherent emission in the current hemisphere from an arbitrary number of soft, energy 
ordered wide-angle gluons in the target hemisphere. 

Additionally, since for the non-photon axis variables one has to resum soft /collinear gluons 
off a hard configuration involving a wide-angle (wrt 7* axis) current jet, we must consider 
a situation where there is a lone gluon at any angle in the current region which also 
corresponds to y = 0. Then one has to resum gluon emissions of this gluonic jet to NLL 
accuracy. 

4. Resummed results and form factors 

For the jet-mass distribution and event shape cross section one obtains an identical formula 
to that derived in Q. For the shape cross section defined as 



""^It has recently been realised that the thrust wrt photon axis and normalised to the total current 
hemisphere energy is a non-global observable in the sense that the emissions in the target and current 
hemisphere are not on the same footing. 
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we have 

R{p) = (1 + asCl ® q + asC^ 5)5(asi)Sg(a„ L) + a,Cf q ^^(a,, L) , (4.2) 

where L = In^. The form factor Sg, which represents resummation off a hard quark 
projectile in Tic, is required tih single logarithmic accuracy and can be obtained from [l| 

Sg(a„L) = / Jq{k')dk' (4.3) 
Jo 

and the function Jq{k'^) has been obtained using the coherent branching technique Q. 
Similarly we get the form factor T,g which represents resummation of a hard gluon jet in 
Tic- S represents the wide-angle soft gluon resummation detailed in This function has 
not been calculated analytically so far but a parametrisation of it in the large Nc limit 
was obtained by fitting to Monte-Carlo results Hence one may expect corrections to 
this piece which are O ^j^r^ (10 % level). The various constants indicate contributions 
from incoming quarks or gluons as indicated by the suffix a and convolutions with the 
appropriate structure functions (pdfs q and g). The upper index denotes the hard parton 
in Tic off which one resums soft emissions. 

For the thrust and the C parameter resummed results can be derived by noting that 
they are proportional to the jet mass in the region of interest and one can use the effective 
relations C = 12p = 6(1 — T) to NLL accuracy. 

Before these results can be compared to experimental data from HERA one has to 
match the results with NLO estimates to get the best possible description over the entire 
range of values of the variables. This procedure is described extensively in ||^]. 

In order to account for hadronisation effects the matched distributions are then shifted 
hy V ^ V — ayV, applying the concept of power corrections |^]. The coefficients ay 
are calculable, while "P is a universal function proportional to 1/Q which depends on a 
non-perturbative parameter oq as well as on Os- The quantity ayV is exactly the same 
hadronisation contribution which appears in the corresponding mean values In fact the 
analysis of the HI data § shows that the Q dependence of the above defined event shape 
means can be consistently fitted by a common parameter ckq — 0.5 and reasonable values 
of as{Mz)- However, fits to event shape spectra using fixed order NLO calculations give in 
general inconsistent parameters (larger spread) compared to those obtained from the means 
Q. One reason for the discrepancies may be missing higher order QCD calculations. The 
hope is that this situation will improve with the present resummations. 

Preliminary comparisons with HI data Q on thrust r = 1 — T and the jet mass p 
spectra over a large range of Q values are shown in figure |l[ The resummed distributions 
are supplemented by power corrections assuming as{Mz) = 0.118 and oq = 0.50. The 
data can be reasonably well described over large parts of the spectrum, the regions of 
agreement extend with rising energy Q. The improvements of the resummations compared 
to fixed order calculations are mainly at small values of the event shape variables. The 
description deteriorates considerably for lower Q data. This can be understood in terms of 
the increasing importance of subleading corrections and huge hadronisation effects. 
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H1 DATA + PRELIMINARY THEORY 
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Figure 1: Differential distributions of thrust t = 1 — T (left) and jet mass p (right). HI data ||^ are 
compared to resummed QCD calculations including power corrections assuming as[Mz) — 0.118 
and So = 0.50. The spectra given at {Q) = 7.5 GeV, 8.7 GeV, 15 GeV, 17.8 GeV, 23.6 GeV, 
36.7 GeV, 57.7 GeV and 81.3 GeV (bottom to top) are multipUed by factors of 10" (n = 0, . . . , 7) 

Summary 

Resummed results are now available for all commonly studied event shape spectra in DIS. 
Preliminary comparisons with data after carrying out matching to NLO and adding power 
corections look very encouraging. A systematic analysis of several DIS event shape distri- 
butions including the dependence on scales, matching schemes, etc. is under preparation. 
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